1 Overview

This short document provides some details on the porpoise movement model developed as an alternative to assess harbour porpoise exposure to wind farms in the North Sea.

Code and data to reproduce the analyses summarized in this document can be found in this repository.

2 Description of approach

Simplified the approach presented in Joy et al. (2022).
We simulate the movement of n_agents from March 1 to Nov 30, following the whale densities.
In general terms, there are four steps of the approach:

  • Step 1: for each agent, choose an initial position in day t.

  • Step 2: given position in day t, propose candidate positions for day t+1 by using random draws of step length and turning angle

  • Step 3: calculate weight for each member of the forecast ensemble. Calculate the probability of each candidate position for each agent

  • Step 4: sample from the set of candidate positions, with the probability in step 3.

2.1 Whale density maps

Density maps are from Gilles et al. (2016). Anita Gilles provided the shape files to Magda Chudzinska. They present seasonal density maps of harbour porpoise in the Noth Sea. Note that they cover spring (March-May), summer (Jun-Aug), fall (Sep-Nov). Therefore our movement model is restricted to March-November.

2.2 Initial positions

To make sure we have porpoises distributed across the study area in the first step, I divided the space into 10 Voronoi blocks of equal area, calculated the expected proportional abundance in each of the blocks and allocated the agents to each of the blocks proportionally to the expected abundance. Within the blocks, I randomly chose the location of the agents, with the probability of the location proportional to the density of each pixel.

2.3 Candidate positions

2.3.1 Step length and turning angles

Nabe-Nielsen et al. (2013) presented distributions of step lengths and turning angles at a 30 minute resolution. The original data was collected at a frequency of 6 seconds in Wilson et al. (2007). Without having the original dataset, I am not sure how to scale up the distribution presented at 30-minute interval to a daily time scale.

I did a quick search and found three tagging studies for harbour porpoise:

  • Read and Westgate (1997): tagged 9 porpoises in the Bay of Fundy and Gulf of Maine - data is available for 8 porpoises

  • Nielsen et al. (2013): tagged two female porpoises in West Greenland - data available for both

  • Sveegaard et al. (2011): tagged 64 harbour porpoises between 1997 and 2007, in order to determine the distribution in the North Sea, the western Baltic, and the waters in between. This dataset is used in several publications, but I could not find the data itself.

I harvested the data from the two papers that have available data, and compare below their distributions to the distribution of step length and turing angles presented in Stalder et al. (2020) (subset of 57 porpoises from the Sveegaard et al. (2011) study).

2.3.1.1 Stalder et al. 2020 data

Note that the data is data is:

  • at a daily time step
  • discriminated by behaviour:
    • Area-restricted movement
    • Directed movement
    • Uncertain

we will not discriminate by behaviour in this exercise, and therefore we only need one distribution that should resemble the mean daily movement

2.3.1.2 Harvested data

I created these distributions using data from Read and Westgate (1997) and Nielsen et al. (2013).

2.3.1.2.1 Step length

I extracted the mean and var of daily displacements from these data, and calculated the parameters of a gamma distribution as:

\[scale = Var(x)/E[x]\] \[shape = E[x]^2/Var(x)\] Using these parameters, I drew random step lengths from a gamma distribution.

2.3.1.2.2 Turning Angle

I drew turning angles from a von Mises distribution. Given that the time step is daily, I allowed the agents to go in any direction (360°).

2.3.1.3 Comparison

At first sight, it seems like the distribution created from the harvested data resembles the distribution presented in Stalder et al. (2020). However, please note this passage from Nielsen et al. (2013): “The average daily travel rates for the two harbour porpoises were 24 km (range 0.6 – 103 km, SD=21.9) and 42 km (range 1.9 – 118.7 km, SD=25.2), respectively, which are in the same range of harbour porpoises from Bay of Fundy (range 13.9 – 58.5 km, SD=13.9), Canada (Read & Westgate 1997) but slightly higher NAMMCO SC/20/HP/08 Harbour porpoise working group November 20133 than found for the North Sea (20.4 kmd-1) and the Baltic Sea (15.4 km d-1) (Sveegaard et al. 2011).”

2.3.1.4 Conclusion

I used the distribution from the harvested data from Read and Westgate (1997) and Nielsen et al. (2013), the “combined” category in the plot above.

2.3.2 Spatial constraints

Candidate positions were considered valid only if:

  • interpolated density at the position was > 0, and
  • was within the North Sea - I am clipping the points using the seasonal density maps as masks

2.4 Calculate weight for each member of the forecast ensemble

For each candidate position, I obtained (using linear interpolation) the expected whale density at the point. The weight for each candidate position was then calculated as: \(weight_{i,t} = density_{i,t} / density_{t-1}\), i.e. the weight of each candidate i at time t is the whale density at that position, divided by the whale density at the accepted position in day t-1.
For each each candidate position, obtain the probability that it will be accepted, by weighting the \(weight_{i,t}\) by the sum of weights:
\(p_{i,t} = weight_{i,t} / \sum_{i=1}^{n.candidates} weight_{i,t}\)

2.5 Choose candidate

For each day and agent, draw one candidate position from the pool of candidate positions, with probability \(p_{i,t}\).

2.6 Practical considerations

I started the agents on March 1, and allowed them to move until Nov 30 at daily steps. For the first day of each season, I changed the underlying density map that the agents followed.

2.7 Parameters

n_agents: 1000
n_candidates: 100

3 Output

3.1 Daily distances travelled

this plot shows the distances travelled by the agents

3.2 Seasonal comparisons of input densities and simulated agent densities

3.2.1 Spring

3.2.2 Summer

3.2.3 Fall

3.3 Animated tracks

Each of the following plots shows the tracks of 4 agents chosen at random from March 1 to November 30. The label “Date” on top of each plot refers to the day of the year, each animation starts on Date 60 (March 1), and finishes on Date 334 (Nov 30).

# References

Gilles, A., S. Viquerat, E. A. Becker, K. A. Forney, S. C. V. Geelhoed, J. Haelters, J. Nabe-Nielsen, et al. 2016. “Seasonal Habitat-Based Density Models for a Marine Top Predator, the Harbor Porpoise, in a Dynamic Environment.” Ecosphere 7 (6): e01367. https://doi.org/https://doi.org/10.1002/ecs2.1367.
Joy, Ruth, Robert S. Schick, Michael Dowd, Tetyana Margolina, John E. Joseph, and Len Thomas. 2022. “A Fine-Scale Marine Mammal Movement Model for Assessing Long-Term Aggregate Noise Exposure.” Ecological Modelling 464: 109798. https://doi.org/https://doi.org/10.1016/j.ecolmodel.2021.109798.
Nabe-Nielsen, Jacob, Jakob Tougaard, Jonas Teilmann, Klaus Lucke, and Mads C. Forchhammer. 2013. “How a Simple Adaptive Foraging Strategy Can Lead to Emergent Home Ranges and Increased Food Intake.” Oikos 122 (9): 1307–16. https://doi.org/https://doi.org/10.1111/j.1600-0706.2013.00069.x.
Nielsen, Nynne .H., Rikke G. Hansen, Jonas Teilmann, and Mads Peter Heide-Jørgensen. 2013. “Extensive Offshore Movements of Harbour Porpoises (Phocoena Phocoena).” Harbour porpoise working group NAMMCO SC/20/HP/08.
Read, A. J., and A. J. Westgate. 1997. “Monitoring the Movements of Harbour Porpoises (Phocoena Phocoena) with Satellite Telemetry.” Marine Biology 130 (2): 315–22. https://doi.org/10.1007/s002270050251.
Stalder, D., van Beest FM, S. Sveegaard, R. Dietz, J. Teilmann, and J. Nabe-Nielsen. 2020. “Influence of Environmental Variability on Harbour Porpoise Movement.” Mar Ecol Prog Ser 648: 207–19. https://www.int-res.com/abstracts/meps/v648/p207-219/.
Sveegaard, Signe, Jonas Teilmann, Jakob Tougaard, Rune Dietz, Kim N. Mouritsen, Geneviève Desportes, and Ursula Siebert. 2011. “High-Density Areas for Harbor Porpoises (Phocoena Phocoena) Identified by Satellite Tracking.” Marine Mammal Science 27 (1): 230–46. https://doi.org/https://doi.org/10.1111/j.1748-7692.2010.00379.x.
Wilson, Rory P., Nikolai Liebsch, Ian M. Davies, Flavio Quintana, Henri Weimerskirch, Sandra Storch, Klaus Lucke, et al. 2007. “All at Sea with Animal Tracks; Methodological and Analytical Solutions for the Resolution of Movement.” Deep Sea Research Part II: Topical Studies in Oceanography 54 (3): 193–210. https://doi.org/https://doi.org/10.1016/j.dsr2.2006.11.017.